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Abstract 

We report a significant improvement of an approximate method of includ- 
ing electron Coulomb distortion in electron induced reactions at momentum 
transfers greater than the inverse of the size of the target nucleus. In partic- 
ular, we have found a new parametrization for the elastic electron scattering 
phase shifts that works well at all electron energies greater than 300 MeV. As 
an illustration, we apply the improved approximation to the (e, e'p) reaction 
from medium and heavy nuclei. We use a relativistic "single particle" model 
for (e, e'p) as as applied to '^^^Pb{e, e'p) and to recently measured data at CE- 
BAF on ^^0(e, e'p) to investigate Coulomb distortion effects while examining 

the physics of the reaction. 
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Electron scattering has long been acknowledged as a useful tool for investigating nuclear 
structure and nuclear properties, especially in the quasielastic region. One of the primary 
attributes of electron scattering as usually presented is the fact that in the electron plane- 
wave Born approximation, the cross section can be written as a sum of terms each with a 
characteristic dependence on electron kinematics and containing various bi-linear products 
of the Fourier transform of charge and current matrix elements. That is, various structure 
functions for the process can be extracted from the measured data by so-called Rosenbluth 
separation methods. The trouble with this picture is that when the Coulomb distortion 
of the electron wavefunctions arising from the static Coulomb field of the target nucleus is 
included exactly by partial wave methods, the structure functions can no longer be extracted 
from the cross section, even in principle. 

Electron Coulomb distortion in elastic and inelastic scattering for various processes has 
been included with various approximations in the past In the early 90's Coulomb 

distortion for the reactions (e, e') and (e, e'p) in quasielastic kinematics was treated exactly 
by the Ohio University group p|^TIl[ using partial wave wave expansions of the electron 
wave functions. Such partial wave treatments are referred to as the distorted wave Born 
approximation (DWBA) since the static Coulomb distortion is included exactly by numeri- 
cally solving the radial Dirac equation containing the Coulomb potential for a finite nuclear 
charge distribution to obtain the distorted electron wave functions. The induced transi- 
tion by a virtual photon is included to first order (the Born Approximation). While this 
calculation permits the comparison of various nuclear models to measured cross sections 
and provides an invaluable check on various approximate techniques of including Coulomb 
distortion effects, it is numerically challenging and computation time increases rapidly with 
higher incident electron energy. And, as noted above, it is not possible to separate the cross 
section into various terms containing the structure functions and develop insights into the 
role of various terms in the charge and current distributions. 

In all of our DWBA investigations of (e, e') and (e, e'p) reactions in the quasielastic 
region, we used a relativistic treatment based on the a — uo model for the nucleons involved. 
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In particular, for the (e, e'p) reaction we use a relativistic Hartree single particle model for a 
bound state [jll| and a relativistic optical model for an outgoing proton |[T^ combined with 
the free space relativistic current operator 

J'' = r + ^^^'"du. (1) 

Using this model, we compared our DWBA calculations with experimental data measured at 
various laboratories for (e, e') and for (e, e'p) jsHTOl and have found excellent agreement 



with the data. We concluded that the relativistic nuclear models are in excellent agreement 
with the measured data and do not need to invoke meson exchange effects and other two- 
body terms in the current that are necessary in a Schrodinger description that uses a non- 
relativistic reduction of the free current operator . Therefore, in this brief report we will 



continue to use our relativistic "single-particle" model to investigate Coulomb distortion 
effects and to compare to the newly measured data from CEBAF. 

To avoid the numerical difficulties associated with DWBA analyses at higher electron 
energies and to look for a way to still define structure functions, our group [0,0] devel- 



oped an approximate treatment of the Coulomb distortion based on the work of Knoll 
and the work of Lenz and Rosenfelder [0. Knoll examined approximations of the Green 
function valid for large momentum transfers (that is, valid for qR > 1 where R is the size of 
the target) while Lenz and Rosenfelder constructed plane-wave-like electron wavefunctions 
which included Coulomb distortion effects. We were able to greatly improve some previ- 



ous attempts along this line |[18| , p!9[| where various additional approximations were made 
which turned out not to be valid. We did have the advantage of having the exact DWBA 
calculation available for incident electron energies up to 400 — 500MeV for checking our 
approximations. We compared our approximate treatment of Coulomb distortion to the 
exact DWBA results for the reaction (e, e'p) and found good agreement (at about the 1-2% 
level) near the peaks of cross sections even for heavy nuclei such as "^^^Pb. The agreement 
was not so good away from the peaks. 

As discussed in our previous papers 0,0], one of the ingredients of our approximate 



electron wavefunction is a parameterization of the elastic scattering phase shifts in terms of 
the angular momentum. In this paper, we briefly review our previous approximation of the 
Coulomb distorted electron wavefunction and present a greatly improved parametrization of 
the phase shifts which works well at all incident electron energies greater than 300 MeV. In 
addition, we will compare our relativistic "single-particle" model to new experimental data 
from CEBAF. 

Our approximate method of including the static Coulomb distortion in the electron 
wavefunctions is to write the wave functions in a plane- wave-like form [ll5| ; 



vl;±(r) = ^ e±^^(L') e*^ e'^'^->' u, , (2) 
P 

where the phase factor 5(L^) is a function of the square of the orbital angular momentum 
operator, Up denotes the Dirac spinor, and the local effective momentum p'(r) is given in 
terms of the Coulomb potential of the target nucleus by 

p'(r) = (^p- ^^V(r)rfrj p . (3) 

We refer to this r-dependent momentum as the Local Electron Momentum Approximation 
(LEMA). The ad — hoc term A = a[p'(r)-r]L^ denotes a small higher order correction to the 
electron wave number which we have written in terms of the parameter a = —aZ{ ^^^^^^^ Y . 
The value of 16 MeV/c was determined by comparison with the exact radial wave functions 
in a partial wave expansion. 

The elastic scattering phase shifts are labelled by the Dirac quantum number k which 
takes on plus minus integer values beginning with one. The eigenvalues of are j{j + 1) 
which equals n"^ — The basic idea of our approximation is to calculate the elastic scattering 
phases and fit them to function of k^. Then to replace the discrete values of /t^ with 
the total angular momentum operator which we then replace by the orbital angular 
momemtum operator since the low k terms where the difference between j and I is 
significant contribute very little to the cross section. The removal of any spin dependence 
apart from what is in the Dirac spinor Up is crucial for defining modified structure functions. 



Based on earlier work by others we fitted the elastic scattering phases shifts to a power 
series in up to second order; 

= fco + &2«^ + &4«^ , (4) 

where the coefficients, 60 ; ^2, ^4 are extracted from a best fit for k values up to about 3pi? 
where R is the nuclear radius. Note that this procedure requires calculating the elastic 
scattering phase shifts for the incident and outgoing electron energies up to k values of 
order 3pi?, which for high electron energies can be quite demanding computationally. We 
refer to these phases as the «;^-dependence phases. This fit to the phases worked very well 
for K values up to approximately k — 3pit!~35 at medium or low energy, but did not fit the 
exact phases shifts very well for higher energies where k — 3pR>50 or more. Since we were 
primarily looking at electron energies in the 300-600 MeV range in our previous work, this 
discrepancy did not present a significant problem. 

However, with CEBAF type energies we need a fit to the phases that will work at any 
incident energy where the overall approximation can be used; that is, for incident electron 
energies greater than about 300 MeV and processes with momentum transfers greater than 
about 1/R. In addition, we would like to avoid calculating all of the elastic phase shifts, 
particularly the very high ones. A reasonable solution is to make use of the fact that the 
higher k phase shifts appoach the point Coulomb phases which have a simple analytical 
form at high energy. At the other extreme, the low k phases corresponding to orbitals which 
penetrate the nucleus are linear in which was the basis or our initial parametrization. 
The difficult phases to fit correspond to k values of order pR which, from a classical point 
of view, corresponds to scattering from the nuclear surface. Moreover, it is well known 
that in electron induced reactions the spatial region around the surface gives the largest 
contribution to the cross section, so it is important to fit the intermediate range as well as 
possible. 

Another goal is to reduce the computer time needed, so we decided to seek a parametriza- 
tion of the elastic scattering phases shifts in terms of which has the correct large 
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behaviour and becomes linear in at low angular momentum. Since we have the correct 
large k behaviour, we need only calculate the exact scattering phase shifts for k values of up 
to of order pr. The large k and small k behaviour are quite different, so we chose to write 
the expression for the phase shift as the sum of two terms with an exponential factor which 
suppresses one of the terms at small k, values and the other at large k values. After some 
experimentation, we find that the following parametrization of elastic scattering phase shift 
describes the exact phase shifts very well: 

2 2 

S{K) = [ao + a2j^]e-fM^ 
aZ k2 

— ^(1 - e (^)xln(l + K^) , (5) 

where p is the electron momentum and we take the nuclear radius to be given hj R = 
1.2^4-'^/'^ — We fit the two constants Oq and 02 to two of the elastic scattering 

phase shifts (k = 1 and k = Int{pR) + 5). To a very good approximation, oq = 5(1) and 
a2 = 4:6{Int{pR) + 5) + aZln{2pR), where Int{pR) replaces pR by the integer just less 
than pR. Note that this parametrization only requires the value of the exact scattering 
phase shift for k = 1 and k = Int{pR) + 5. As shown in Fig. |I|, the K^-dependence phase 
parametrization breaks down for high k values and has large deviations for mid-range k 
values. The new phase parametrization fits the exact phases very well for electron energy 
of E = 2400 MeV on ^^O, although the new phase parametrization does still show some 
small deviations from the exact phases for n values around 20 to 30 which is in the surface 
region. Clearly additional terms could be added to the parametrization to obatin a better 
fit. However, as we shall see below, the simple fit that we have used reproduces the cross 
section quite well. 

Using the new phase shift parametrization and the local effective momentum approxima- 
tion we construct plane-wave-like wave functions for the incoming and outgoing electrons. 
Since the only spinor dependence is in the Dirac spinor all of the Dirac algebra goes through 
as usual and we end up with a M0ller-like potential which contains an r-dependent momen- 
tum transfer. It is then straightforward to calculate the (e, e'p) cross sections and modified 
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structure functions. Please see our previous papers |ll^ , [T5[] for details. 

In most (e, e'p) experiments, there is sufficient energy resolution that protons knocked 
out of different shells can be examined. It is common to report the experimental results 
in terms of the reduced cross section pm as s function of missing momentum prn-, which is 
defined by 

, , 1 d?a . , 

P"^^""^ ^ PEpa.p dEfd^fdnp ' ^ ^ 

where the missing momentum is determined by the kinematics pm = P — q where P is the 
outgoing proton momentum and q is the asymptotic momentum transfer from the electron 
defined by q = Pi — p/. For plane wave protons in the final state pm is related to the 
probability that a bound proton from a given shell have momentum p^. For the off-shell 



electron-proton cross section a^p we use the form 'ccl' given by de Forest pO[. For distorted 
outgoing protons, this reduced cross section is just a convenient way of comparing experiment 
and theory since the theory results for the cross section can have the same factors removed. 
Note that all calculations will be carried out in the laboratory system (target fixed frame). 

While there are two experimental kinematic arrangements commonly used in (e, e'p) 
experiments with designations of parallel kinematics and perpendicular kinematics, in the 
present work, we consider only perpendicular kinematics. In perpendicular kinematics, 
the momentum transfer q is held fixed along with the magnitude of the momentum of 
the outgoing proton while the angle between q and P is varied. The calculated reduced 
cross section is compared (by means of a linear least squares fit) to the similarly reduced 
experimental cross section to extract an overall scale factor which is the spectroscopic factor. 
The spectroscopic factor contains two factors, the occupation probability of a proton in a 
given orbit and the overlap of the residual nucleus with the A — 1 nucleons in the target. 

As a test case, we calculate the reduced cross sections with the new phases for a heavy 
nucleus, ^°^P6. Figure ^ shows the reduced cross section as a function of the missing mo- 
mentum Pm for knocking protons from the 3si/2 orbital of "^^^Pb. The incident electron 
energy Ei = 412 MeV, and the outgoing proton kinetic energy is Tp = 100 MeV. We have 



chosen P = q which corresponds to an electron scattering angle of 9e = 74°. The solid 
line is the result of the full DWBA [Q, the dashed curve is the result with the new phase 
shift parametrization, and the dotted curve is the result with K^-dependence phase shift 
parametrization. The dashed curve obtained by using the new phases clearly reproduces the 
exact result much better than the previous K^-dependence phase parametrization over the 
whole region. 

We also apply the new phase shift paramtetrization for the case of high energy electron 
scattering on the light nucleus ^^O where protons are removed from the pi/2 and J93/2 orbits. 
The incident electron energy Ei = 2441.6 MeV and the outgoing proton kinetic energy 
Tp = 427 MeV as shown Fig. ^. In this figure, the solid curves are the approximate DWBA 
results using the new phase shift parametrization, the dotted curves are the PWBA results 
without Coulomb distortion, and the data are newly measured from CEBAF as reported in 



the dissertation of Gao [^. Note that our exact DWBA code cannot evaluate such high 
energy processes without extensive modification which we have not done. 

As expected, the effect of Coulomb distortion on such a high energy electron induced 
process is very small except possibly at large missing momentum. Note that the Coulomb 
effects for ^^O in the medium energy region (500 MeV) was of the order of 3% |TD[. This 



fit to the experimental data using our relativistic "single particle" model for the nucleon 
wavefunctions results in spectroscopic factors of 61% for the pi/2 orbital and 70% for the 



P3/2 orbital. In our analyis of Saclay data p3[ at lower electron energies using a similar 



nuclear model we found spectroscopic factors of 54% and 57% respectively [|T0|- 

In summary, we have improved our previous approximate method of including Coulomb 
distortion effecs in (e, e'p) reactions from nuclei. The improvement involves a better 
parametrization of the elastic scattering phases shifts which has the correct behaviour for 
large angular momenta and requires the calcuation of only two phase shifts (for k = 1, 
and for k equal to Int{pR) + 5). We showed that even for (e, e'p) on "^^^Pb the cross sec- 
tion calculated with our approximation using the improved parametrization of the phase 
shifts agrees with the exact DWBA result quite well even out beyond the second maxima. 
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This is a significant improvement over our previous approximation for the phase shifts. In 
addition, we compared our relativistic "single-particle" model for (e, e'p) from ^^O to the 
recently measured cross section at Thomas Jefferson Lab and found excellent agreement for 
the removal of a proton from the pz/2 and pi/2 shells with reasonable spectroscopic factors. 

Our improved approximate method of including Coulomb distortion in electron scattering 
reactions works for high energy electrons as well as for more moderate energies (300 — 
SOOMey), and for experiments at the few percent level this approximate way of including 
Coulomb distortion is adequate. More importantly, as discussed in our previous paper, this 
"plane-wave-like" approximation permits the extraction of "structure functions" even in the 
presence of strong Coulomb effects and thus provides a very good tool for looking into the 
response of the nucleus to "longitudinal" and "transverse" photons. 
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FIG. 1. The comparison between the exact, K^-dependent, and new phases in ^^O for 
i^max = 100 and energy E = 2441 MeV. The diamonds are the exact phases, the dashed curve for 
fit, and the soUd curve for the new phase shift parametrization. 
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FIG. 2. The reduced cross section for ^^^P6 for the 3si/2 orbit with perpendicular kinematics. 
The kinematics are Ei = 412 MeV, and proton kinetic energy Tp = 100 MeV. The sohd hne 
is the full DWBA result, the dashed line is the approximate DWBA using the new phase shift 
parametrization, and the dotted line is the approximate DWBA with the K^-dependent phases. 
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FIG. 3. The cross sections for O from pi/2 and orbits for perpendicular kinematics. 
The kinematics are Ei = 2441.6 MeV, proton kinetic energy Tp = 427 MeV, and energy transfer 
Lv = 436 MeV. The sohd hnes are the new approximate DWBA results, the dotted lines are the 
PWBA results, and the data are from CEBAF [||]. 
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